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ABSTRACT 

t-H , 

I Spatial configurations of randomly charged polymers, known as polyampholytes 

. (PAs), are very sensitive to the overall-jxcess charge Q. Analytical arguments, 

| supported by Monte Carlo simulationsErQ and exact enumeration studies,H lead to 

00 ■ the following picture: For Q < Q c w qoVN (qq is the elementary charge, TV is the 

C> number of monomers in the polymer), the radius of gyration R g of the polymer 

decreases with decreasing temperature T and the polymer becomes compact, 
while for Q > Q c the polymer stretches with decreasing T. At low T, the dense 
states are described by Debye-Hiickel theory, while the expanded states resemble 
a necklace of globules connected by strings. At such temperatures, the transition 
between the dense and the expanded states with increasing Q, is reminiscent of 
the breakup of a charged drop. 

O - 

O ■ 1. Introduction 



Statistical mechanics of randomly charged polymers, called polyampholytes (PAs), 
is a challenging subject because it embodies an interesting combination of long range 
interactions and randomness. While the physics of homogeneous polymers has a rea- 
sonably firm basis,0 considerably less is known about heteropolymers, although the 
latter present an extremely rich problem of biological significance]! In this presenta- 
tion we review some properties of PAs. 

The spatial extent of a polymer is characterized by the critical exponent is, which 
relates its radius of gyration (r.m.s. size) R g to the number of monomers N, by the 
power law R g oc N u . The polymer will be called "compact" if v — 1/d, where d 
is the dimension of the embedding space, and "stretched" if is = 1. The simplest 
model of PAs is as a flexible chain of monomers, each of which has a charge ±go 
selected from a well defined ensemble of quenches. The polymer has a characteristic 
microscopic length scale a, such as the range of excluded volume interactions, or the 
nearest neighbor distance along the chain. The monomers of the PA interact both 
via (short range) excluded volume interactions and long range (unscreened) Coulomb 



interactions. In the simplest ensemble of quenches, each monomer takes a charge 
<7i = ±g independently of all the others; i.e. q^q] = f^g^, where the overline indicates 
averaging over quenches. While the average excess charge Q = J2i Qi of such PAs is 
zero, a "typical" sequence has Q of about ±Q C , with Q c = qoN 1 ^ 2 . This statement, as 
well as the definition of Q c , is unrelated to the embedding dimension d. However, the 
importance of charge fluctuations (both for the overall polymer, or for large segments 
of it) does depend on the space dimension. The electrostatic energy of the excess 
charge, spread over the characteristic size of an ideal polymer (R g oc iV 1 / 2 ), grows 
as Q 2 /Rg~ 2 ~ jV"( 4-d )/ 2 . This simple dimensional argument shows that for d > 4 
weak electrostatic interactions are irrelevant. (The excluded volume effects are also 
irrelevant in d > 4.) Thus, in d > 4 and at high temperatures, the PA behaves 
as an ideal polymer with an entropy-dominated free energy of the order of —NT. 
However, on lowering temperature it collapses into a dense state, taking advantage 
of a condensation energy of the order of — iVg 2 /a d_2 . 

We are primarily interested in the behavior of PAs in d — 3. We present both 
numerical and analytical evidence that the behavior of PA in d < 4 is extremely 
sensitive to the value of Q. In particular, we show that for Q < Q c the PA assumes 
compact configurations, while for Q > Q c it is in an expanded state. We suggest a 
"necklace" model which qualitatively describes the latter state. 

2. Properties of Polyampholytes 



2.1. Neutral Heteropolymers: Short-Range vs. Long-Range Interactions 

It is interesting to compare and contrast the behaviors of heteropolymers with 
short and long range interactions. Consider a polymer represented by a self-avoiding 
chain carrying a quenched sequence of random charges which interact only at short dis- 
tances, with similar charges repelling each other and the opposite charges attracting.! 
In d = 3, at high temperatures T, the presence of such short range interactions does 
not influence the behavior of the polymer and it behaves as a self-avoiding walk 
(SAW). If Q is small, i.e. the positive and negative charges are almost balanced, 
upon decrease of temperature, charges of opposite sign are more likely to be in the 
vicinity of each other. This introduces an effective attraction which increases in 
strength with decreasing temperature, and at a certain temperature the polymer col- 
lapses into a globule, as in a regular ^-transition. When the imbalance Q between the 
two types of charges increases, the transition temperature decreases, as demonstrated 
in Fig. ^|a, and for sufficiently large imbalance the transition disappears altogether.! 
Since the behavior of the system varies very gradually with increasing Q/N, it does 
not really matter whether the ensemble of chains consists only of quenches with ex- 
actly vanishing Q, or is obtained by independently choosing each qi = ±q with equal 
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Fig. 1. Phase diagrams for random polymers with (a) short-range, and (b) Coulomb interactions, 
in d — 3. Note the different scaling of the excess charge Q in the two figures. 



probability. In the latter case Q may deviate from zero by an amount of order qoN 1 ' 2 ; 
such fluctuations have a negligible contribution to Q/N for large N. 

The long range Coulomb interactions in random PAs are always relevant at d < 4. 
An approximate treatment using a Debye-Hiickel type theory! leads to the conclusion 
that a PA minimizes its free energy by assuming a configuration in which positive 
charges are predominantly surrounded by negative charges, causing screening as in 
a regular electrolyte. A detailed calculation suggests!'! the following picture: At a 
finite T, the PA can be divided into segments containing n monomers each, where 
n defines the length scale at which electrostatic interactions become of order T, i.e. 
Q 2 (n)/R g (n) ~ qQn/(an u °) = T in d = 3, where u w 0.59. The overall configuration 
is then composed of n-mer blobs: Inside the blob the PA behaves as an uncharged 
polymer, while the blobs form a compact structure resembling an electrolyte. The 
resulting PA is compact (i.e., has R g ~ iV 1//3 ) at all temperatures, with a T-dependent 
density. On the other hand, scaling arguments^ assume that the ^-dependence of the 
R g of a random PA is determined from the relation Q 2 /R g ~ q^N/Rg w T, leading 
to the conclusion that R g ~ iV in d = 3. This apparent contradiction is resolved by 
notingi that the spatial configurations of PAs are extremely sensitive to the overall 
charge Q. In particular, it can be shownllfl that there is a critical charge Q c = qoN, 
such that for Q < Q c the PAs are compact, while for Q > Q c they are stretched. 
The diagram in Fig. |l]b indicates such separation between contracted and expanded 
regimes of the PA. In the following subsections we shall present different arguments 



(a) (b) 
Fig. 2. Low T configuration of 64-monomer PAs with (a) Q = 12, and (b) Q = 0. Dark and bright 
shades indicate opposite charges. 



and numerical results supporting this conclusion. 

2.2. Poly ampholytes with Excess Charge: High Temperature Arguments 

For d < 4, electrostatic interactions are relevant and the high temperature phase 
is no longer a regular self-avoiding walk. At high T, the behavior of the polymer 
can be explored perturbatively. For the ensemble of uncorrelated charges (<Mj = 
0), the lowest order (1/T) correction to the quench-averaged R 2 vanishes.^ How- 
ever, restricting the ensemble to yield fixed Q, slightly modifies the quench-averaged 
charge-charge correlations. In particular, the two-point correlation function becomes 
qjq] = (Q 2 — Q 2 )/N 2 for i ^ j. This small (order of 1/N) correction to the correlation 
function causes a significant change in R 2 due to the long range nature of the Coulomb 
interaction. A 1/T-expansionll has a first order term proportional to Q — Q c . Thus 
the size of a PA tends to decrease upon lowering temperature if Q is less than Q c , 
and increases otherwise. 

Confirming the 1/T-expansion, a strongQ-dependence of R g has been seen in 
Monte Carlo simulations:§ It has been shownl that for Q < Q c the PA contracts with 
decreasing T, while for Q > Q c it expands. Fig. || depicts low temperature equilibrium 
configurations of two PAs, one of them has Q slightly larger than Q c and is strongly 
stretched, while the other has Q = and forms a globule. The two behaviors are 
separated by the vertical line at the top of the phase diagram in Fig. |I|b. It should 
be noted that the arguments leading to location of this line are independent of d, as 
long as d < 4. 



2.3. Poly ampholytes with Excess Charge: Low Temperature Arguments 



Monte Carlo results indicate that PAs with Q = form dense globules at low 
temperatures. We can use this observation as a starting point for the investigation of 
the dependence of low-T shapes on Q in d = 3. The energy (or rather the quench- 
averaged free energy) of the PA is phenomenologically related to its shape by 

bO 2 

E = -e c N + 1 S+°-^-. (1) 

The first term is a condensation energy proportional to the volume (assumed com- 
pact) and e c ~ q 2 /a, the second term is proportional to the surface area S (with a 
surface tension 7 ~ q^/a 3 ), while the third term represents the long-range part of 
the electrostatic energy due to an excess charge Q (b is a dimensionless constant of 
order unity). The correctness of separating the energy into the 3 terms appearing 
in Eq. [1] is not self-evident in random systems. A justification for the assumptions 
implicit in this equation is provided by an exact enumeration study of the ground 
states of PAs defined on a latticel Fig. [I] depicts the ground state energies of all 
2080 possible quenches for 12-step (13-atom) chains. (The horizontal axis represents 
an arbitrary numbering of the quenches.) The energies are clearly separated into 7 
bands, corresponding to excess charges of Q = 1, 3, 5, • • •, 13. (There is only one 
quench with Q = 13.) While each band has a finite width, we see that the energy of 
a PA can be determined rather accurately by only specifying its net charge Q, while 
the additional details of the quench appear only in the width of the band. 

The optimal shape is obtained by minimizing the overall energy in Eq. |. The first 
term is the same for all compact shapes, while the competition between the surface 
and electrostatic energies is controlled by the dimensionless parameter 

= Q 2 = _Q^_ _ Q^_ (9] 

a ~ 167lR 3 7 12V 7 ~ Q\ ' U 

Here, R and V are the radius and volume of a spherical drop of iV particles, and 
we have defined the Rayleigh charge Qr. We note that in the case of PAsii Q R « 
q^N 1 / 2 = Q c . The behavior of a system described by Eq. p] has been analyzed in the 
past in the contexts of charged conducting drops, as well as liquid drop models of 
atomic nuclei. If one restricts the possible shapes of a globule to spheroids, then one 
finds that for sufficiently large values of a (of order one or more), the energy of the 
system can be decreased by distorting the drop into a prolate spheroidal shape. Both 
MC simulations^ and experimental resultsi indicate that for Q < Qr, the R g of PA 
at low temperature is almost independent of Q, while it increases extremely fast as a 
function of Q for larger charges. 
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Fig. 3. Ground state energies of (arbitrarily numbered) quenches. 

The transition between compact and expanded states with increasing Q is rep- 
resented by the vertical line at the bottom of the diagram in Fig. [TJb. The presence 
of this line signifies the instability of the spherical shape, but provides no indication 
of what the stable state of the PA looks like. This question is taken up in the fol- 
lowing section. It should be stressed, that although both high and low temperature 
arguments seem to provide a consistent picture of transition between compact and 
extended states, the equality between Q c and Qr depends on the space dimension. 
For d < 3 (d > 3) Qr increases with iV slower (faster) than Q c , and therefore the 
high-T criterion will not coincide with the low-T limit. 

3. The Necklace Model 

3.1. Uniform Charges 

Although the ellipsoidal globule may have a lower energy than the spherical one, 
the former is not the minimum energy configuration of a charged drop: A uniformly 
charged drop minimizes its energy by splitting into two equal, infinitely removed, 
droplets for a > 0.3, and the number of such droplets increases with a. Obviously, 
the PA must maintain its connectivity. We can constrain the overall shape to remain 
singly connected by linking the droplets via narrow tubes of total length L and diam- 
eter a. As long as La 2 R 3 , most of the charge remains in the spheres. The total 
electrostatic energy is proportional to Q 2 /L, while the surface energy cost grows as 



'jaL; equating the two gives L oc Q. For large Q, the PA will look like a necklace of 
globules connected by narrow strands. The configuration depicted in Fig. indeed 
bears some resemblance to such a necklace. 

It should be noted that R g in this simplified discussion is proportional to the 
typical charge Q, and therefore proportional to N 1 ^ 2 . This serves as an indication that 
the quench-averaged configurations are not compact. The picture has one important 
shortcoming for PAs: it assumes that the excess charge is uniformly distributed along 
the chain and disregards the strong charge fluctuations. It thus applies to such 
problems as an alternating sequence of charge s§ to which some additional charge has 
been added, or to a weakly charged homogeneous polymer in which the attraction 
(and, therefore, the condensation energy) is provided by a different mechanism, such 
as short range attractions. 

3.2. Effects of Randomness 

We now illustrate the difficulties caused by randomness for the case of an unre- 
stricted PA. Since Q 2 = q^N, where the overline denotes an average over the ensemble 
of all quenches, we have a — 1. A uniformly charged drop is unstable to splitting al- 
ready for a ~ 0.3, and thus a typical random PA is expected to form several globules 
connected by narrow tubes. Now consider splitting sequences of N monomers with 
total charge constrained by a particular a into two equal subchains of charges Qi and 
Q 2 . It is easy to show that each segment has a su bcham = (1 + &)/2, while the mean 
product of the charges is Q1Q2 = Q^Ni^a — l)/4. The subchains have, on average, 
values of a close to unity. Also, for a = 1, the average value of the product of charges 
vanishes. We thus have the paradoxical situation in which most spherical shapes are 
unstable, while there is on average no energetic gain in splitting the sphere into two 
equal parts. 

Therefore, charge inhomogeneities drastically modify the necklace picture. The 
resulting PA is probably still composed of rather compact globules connected by a 
(not necessarily linear) network of tubes. The globules are selected preferentially from 
segments of the chain that are approximately neutral (or at least below the instability 
threshold), while the tubes are from subsequences with larger than average excess 
charge. It can be shown0 that the probability of finding a very large neutral segment 
in a random sequence of charges is large, i.e. it is possible to build a configuration 
consisting of one very large neutral globule with highly charged "tails" sticking out 
of it. 

While the necklace model provides a convenient starting picture for the behavior 
of PAs, it does not encompass the entire complexity of the problem. The probability 
distribution of the R g of the ground states of PAs obtained from exact enumeration 
study! has a large peak for small R g , indicating that most of the configurations 
are compact, as well as a slowly decaying tail at large R g . This tail represents the 



expanded configurations, and determines the behavior of the quench-averaged R g . 
Indeed, for iV < 13 we foundi that 80% of all quenches can be classified as compact, 
i.e. their R g increases as A 1 / 3 , while the remaining 20% have quench-averaged R g s 
increasing with TV at least as fast as that of a SAW. 

We can summarize our current knowledge of the behavior of random PAs as fol- 
lows: The exact enumeration and Monte Carlo studies confirm that random PAs have 
(on the average) expanded spatial conformations, although further work is needed to 
verify whether v — 1. The necklace model provides a useful qualitative view of the 
ground state configurations. Further studies are needed to put the model on a more 
quantitative basis. 
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